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Abstract 

The B s —7- Ds*^K^ decays allow a theoretically clean determination of tfi s + 
7, where <fi s is the Bg-B® mixing phase and 7 the usual angle of the unitarity 
triangle. A sizable B s decay width difference Ar s was recently established, which 
leads to subtleties in analyses of the B s — > D^^K* branching ratios but also 
offers new "untagged" observables, which do not require a distinction between 
initially present B® or B® mesons. We clarify these effects and address recent 
measurements of the ratio of the B s — > D^rK^ , B s — > D^rn^ branching ratios. In 
anticipation of future LHCb analyses, we apply the SU (3) flavour symmetry of 
strong interactions to convert the l?-factory data for Bj — > D^^tt^, Bd — > Dfir T 
decays into predictions of the B s — > D^^K^ observables, and discuss strategies 
for the extraction of <p s + 7, with a special focus on untagged observables and the 
resolution of discrete ambiguities. Using our theoretical predictions as a guideline, 
we make simulations to estimate experimental sensitivities, and extrapolate to 
the end of the planned LHCb upgrade. We find that the interplay between the 
untagged observables, which are accessible thanks to the sizable Ar s , and the 
mixing-induced CP asymmetries, which require tagging, will play the key role for 
the experimental determination of <p s + 7. 
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1 Introduction 



The decays B s — > D^^K^ only receive contributions from tree-diagram-like topologies^] 
Since both B® and B° s mesons can decay into the Ds^K^ final states, interference effects 
between B®-B® mixing and decay processes allow a theoretically clean determination of 
the CP- violating phase <p s + 7 [l}[2], where S is the B®-B® mixing phase and 7 the 
corresponding angle of the unitarity triangle. As <p s can be extracted separately, with 
the latest experimental average given by (3) 



7 



s =(-2.5^)°, (1) 

7 can be determined. 

The central question is then whether this value will agree with 7 determinations from 
decays with penguin contributions, such as the B® — > 7r + 7r~, B® — > K + K~ system [I]. 
The current picture of direct determinations of 7 from tree decays can be summarized 
as follows: 

f (66 ± 12)° (CKMfitter Collaboration |5|) 
\ (76 ± 10)° (UTfit Collaboration §). ^ ' 

On the other hand, a recent analysis of the B® — > 7r + 7r~, B° s — > K + K~ system gives 

7 = (68 ± 7)°, (3) 

where the error also takes S"?7(3)-breaking corrections into account [7]. 

In the present paper, we assume that the relevant decay amplitudes are described 

by the Standard Model (SM). Applying the formalism developed in Ref. [2], we shall 

(*)± r 

explore the B s —¥ D s channels both in view of recent experimental developments 

and measurements to be performed by the LHCb collaboration in this decade. 

Using the B Q S — > J/ip<fi channel, the LHCb experiment has recently established a 
non-vanishing decay width difference of the 5 s -meson system, which is characterized by 
the following parameter [8] : 

Vs = 2T7 = L 2Fs " = °' 088 ± °- 014 ' (4) 

Here r^ 1 = T S = [r[ s) + T^j/2 = (0.6580 ± 0.0085) ps" 1 ^ denotes the inverse of the 
B s mean lifetime tb 8 - A discrete ambiguity could also be resolved |9], thereby leaving us 
with the sign of Ar s in Q, which is in agreement with the SM expectation (for a recent 
review, see [To]). 

This new development in the exploration of the S s -meson system has important 
consequences: 

• Untagged B s decay data samples, where no distinction is made between initially, 
i.e. at time t = 0, present B° s or B® mesons, allow for an extraction of interesting 
observables [2]. 



1 We use the notation £> s = 



1 



• A subtle difference arises between the branching ratios extracted experimentally, 
and those usually considered by theory (TTJ. 

First measurements of the B s — > DfK T branching ratios are available from the 
CDF fl2l Belle [i~3l and LHCb flil collaborations: 



RR , R n ±KT\ f 0.097 ±0.018 (stat.) ±0.009 (syst.) [CDF], 

S^'Ar = 0.0651H (stat.) [Belle], (5) 

Dtt^,^w,ir ja, [ 0.0646 ±0.0043 (stat.) ±0.0025 (syst.) [LHCb]; 

the errors of the Belle result are dominated by the small B s — > DfK T data sample. 
We shall clarify the impact of Ar s on this ratio of CP-averaged experimental branching 
ratios and convert the experimental numbers into constraints on the hadronic parameter 
characterizing the interference effects discussed above. 

As was pointed out in Ref. [2), the observables of the B s — > Dl^K^ channels can 
be related to those of the — > D^^ir^ decays through the [/-spin symmetry of strong 
interactions. We shall use 5-factory data for the latter decays obtained by the BaBar 
and Belle collaborations, with further constraints from B^ — > Df7c T modes, to make 
predictions for the B s — > D^^K^ observables that will serve as a guideline for the 
expected experimental picture. In this analysis, we specifically find that - thanks to the 
sizable value of AT S - untagged data samples of B s — )■ D s K T decays can be efficiently 
combined with mixing-induced CP asymmetries of tagged analyses to extract <p s + 7 in 
an unambiguous way. 

The outline is as follows: in Section [2j we discuss untagged measurements of the 
B s — > D^^RT decays and their effective lifetimes, addressing also the results listed in 
(p|. In Section [3j we apply SU(3) flavour symmetry to extract the hadronic parameters 
characterizing the B s — > Di*^ ± K T decays from the 5-factory data for the B c i — > D^ ±r K T 
and Bd — > Dfn T channels. In Section El we discuss the extraction of S + 7 from the 
tagged and untagged B s — >■ D s K T observables, with a special emphasis on resolving 
the discrete ambiguities. The hadronic parameters obtained in Section [3] are used in 
Section [5] to predict the relevant B s — > D^^K^ observables, which then serve as an 
input for exploring the experimental prospects. Finally, we summarize our conclusions 
in Section |6l 



2 Untagged Observables 

The time-dependent, untagged B s — > K~ decay rates can be written as follows [2): 
(T{B a (t) -)• D^ + K-)) = F(B° s (t) -)• D^ + K-) + T{B° s {t) -> D^ + K~) 

= Re- t/TB 



cosh ( y s ] +^. A rsinh (y s — 



(6) 



where 

R = (T(B s (t) -+ DU+K~))\ =n = (T(B s (t) -+ D^K+))\. (7) 



2 



The time-dependent, untagged B s — > K + rate into the CP-conjugate final state 
can be straightforwardly obtained from ^ by replacing Aat with Aat- The latter 
observables take the form 



A Ar = -(-1) 



2x* 



cos 



+ j + 6 a ), Aat 



-v. 



, 2x, 



cos 



+ J-S S 



where L denotes the angular momentum of the final statdj the hadronic parameter 
x s oc Rb quantifies the strength of the interference effects between the B" — >■ D s K~ 
and B° s -> Di* )+ K~ decay processes induced through B®-B® mixing, and <5 S is an asso- 
ciated CP-conserving strong phase difference j2]; the parameter Rt oc | (AV^&) | ~ 0.4 
measures one side of the unitarity triangle. 

The branching ratios of B s decays are determined experimentally as time-integrated 



untagged rates 11 , 15 



BR(B S -+ J DW ± ^) cxp = - / (T(B S -+ D^RT)) dt. 



(9) 



On the other hand, the branching ratio corresponding to the untagged rate at t — 0, 
where B®-B° s mixing is "switched off" , is usually considered by theorists. The conversion 
between this theoretical branching ratio and the experimental branching ratio is given 
as follows [IT] : 



BR(B S -»■ D<f )+ K- 



itheo 



1 + ^Ar y s 



/exp; 



(10) 



where an analogous expression involving A^r holds for the D s *' K + final states. It is 
interesting to note that we have 



/theo 



BR(B S -»■ D^ + K- 
thanks to ([7]), which implies 

BR(B S -> D { * )+ K' 



BR(B S D&-K-* 



/theo 



/cxp 



BR(5 S D ( * ] ~K 



/exp 



1 + ^Ar Us 
1 + -4Ar y s 



(11) 



(12) 



Consequently, an established difference between the experimental B s — > D s K + and 
B s —7- D^ + K~ branching ratios would imply a difference between the ^.Ar and ^.Ar 
observables (see also Ref. [16]): 



BR(5 S ->■ Di* )+ K 



BR{B S -> 



/exp 



/exp 



BR(5 S ->■ J Di* K K-) cxp + BR(5 S ->■ Df>-K+) 



(*)- 



cxp 



-4 



Ar 



af 



2 + y s (^ A r + ^Ar) 



(13) 



In order to relate theory to experiment beyond an accuracy corresponding to the size 
of y s ~ 0.1, we need theoretical input to determine Aaf and A^r- In Section K3l we will 



2 For simplicity, we did not introduce a label to distinguish between DfK and D* + K 
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Figure 1: The colour-allowed tree (a) and exchange (b) topologies contributing to the 
B° s -> D { * )+ K- decay in comparison with the colour-allowed tree (c) topology of the 
577(3)-related B® — > Da 7T~ channel which does not receive exchange contributions 
because of the flavour content of its final state. 



see that this results in large uncertainties for these observables. However, this input can 



be avoided with the help of the effective decay lifetimes 11 , defined as 
f™t(T(B B ^Dl* )+ K-))dt 



r(T(B s ^D { s * )+ K-))dt I-Vi 



l + 2A A rys + yj 
1 + *4. A r y s 



(14) 



with an analogous expression for the lifetimes T e s of the CP-conjugate Ds K + final 
states. We then obtain 



BR(B S D ( * )+ K- 



Jtheo 



[2-(l 



y 2 s )r eS ] BR(B S ^D^ + K- 



lexpj 



(15) 



and correspondingly for the D^~K + final states. These general relations hold also 
should the B s — > Dg^K^ decay amplitudes receive contributions from physics beyond 
the SM, which is not a plausible scenario. 

Let us now have a closer look at the ratio (JsJ) . Since the B® — > D~n + , B® — > Dfn~ 
decays are flavour-specific, their ^4.Ar, -Aat observables vanish. The branching ratios 
entering ^ are averages of the experimental branching ratios over the final states: 



BR(B S -> Df 



exp 



[BR(B S -> DtK-) cxp + BR(B S -+ D;K + ) cxp ] 



(16) 



with an analogous expression for BR(B S — > Dfir T ) cxp . Using and its B s — > D 
counterpart yields 



BR{B S -> D { * ]± KT) 



cxp 



BR{B S Di* ]± n^ 



'cxp 



A Ar + A 



BR{B S ->■ Di* ]± RT) 



theo 



J BR(B S -»■ D^ic* 



(17) 



'theo 



plitudes of the B° s ->■ D^K~ and B° s -> D^tt~ decays 
by experimental data 23 . In Fig. [TJ we illustrate the decay topologies characterizing 



(*)- 



The "factorization" of hadronic matrix elements is ex pect ed to work well for the am- 

22 1, which is also supported 
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these decays. Using the SU(3) flavour symmetry to relate the B° s — > Di*^ + K ampli- 
tude to that of the B° s — > D^ + n~ channel (and correspondingly for the CP-conjugate 
processes), the ratio of the theoretical branching ratios in (17) allows the extraction of 
the hadronic parameter discussed in detail in Ref. |2j: 



.1 u 



'theo 



thco 



1. 



Here 



A 2 



e 



I A- 



0.0534 ± 0.0005 



(18) 



(19) 



involves the Wolfenstein parameter A = |V^ S | 
coefficient can be written in the following form: 



0.2252 ± 0.0009 M, while the CW 



C (*) 



where the $ are straightforwardly calculable phase-space factors, and 



(20) 



/it 



Ik F 



(*) 



(M 



(21) 



describes factorizable 5*?7(3)-breaking corrections through the ratios of decay constants 
f K /f n = 1.197 ± 0.006 (24] and form factors^ On the other hand, the non-factorizable 
SU (3)-breaking corrections affecting the ratio of the colour-allowed tree amplitudes gov- 
erning; the B° s -> D { s * )+ K- and B° s ->■ L>i* )+ vr- channels are described by 



ai(Di* } 7r 



(22) 



Finally, takes into account that the B° s — > D^ + K~ decays receive also contributions 
an 

as can be seen in Fig. [I] 



from exchange topologies, which have no counterparts in the B® — > D\ ; 7T processes 



K~ 



T D (*) 



K- 



+ E D^+K- 



(23) 



Following the phenomenological analysis of Ref. [23] using experimental data to make 
factorization tests and to constrain the exchange topologies, we find M* ] ~ 1-00 ± 0.02 
and A/g* ~ 0.97 ± 0.08. The exchange contributions can be probed further in the 



3 For the calculation of the form- factor ratio in (21) we have assumed that the q 2 dependence is 
identical to that for Bd — > D^*'£u decays 25 
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Figure 2: Compilation of measurements of the ratio of branching ratios as given in ^ 
and comparison with the lower bound in (26). The theoretical prediction indicated by 



the vertical band corresponds to ( 63 ) as given in Section 5 



future through the B° s — > D^ + ti channel, which receives only contributions from such 
topologies pi. Finally, we obtain the numerical value 



C w = 0.67 ±0.05. 



(24) 



Using now mh and (17), we arrive at 



x s = Us cos5 s cos(4> s + 7) 



cm 



BR(B S -> D { * ]± KT 



'cxp 



BR(B S -> D^nT 



'cxp 



l + Us cos 2 S s cos 2 (<p s + 7), (25) 



where x s was defined as a positive parameter pi. For the numerical values of S and 7 in 
g} and Q, respectively, the CDF result in Q gives x s = 0.46 ± 0.27 (BR) ± 0.11 (C) ± 
0.04 (8 S ). This value for x s is consistent with theoretical expectations [2] and the picture 
discussed in the next section. On the other hand, the central values of the LHCb and 
Belle results in (pjl) do not give real solutions for x s . The requirement that the argument 
of the square-root in (25) is positive can be converted into the following lower bound: 



BR{B S -> D { ; ]± RT) 



cxp 



> 



'cxp 



CM 



1 - y 2 s cos 2 5 S cos 2 (4> s + 7) 



0.080 ±0.007, 



(26) 



which is shown in Fig. [2] We observe that the LHCb result for the ratio of branching 
ratios would need to increase by about two standard deviations to satisfy this bound 
and to give a real solution for x s . 

In the next section, we shall use data from the B factories to obtain a sharper picture 
of the hadronic parameters, including the CP-conserving strong phases 5 S . 
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3 Hadronic Parameters from Bd — >• iD/*? =l= 7r =F Data 

a (s) 

Using the [/-spin flavour symmetry of strong interactions, the hadronic parameters x s 
and 5 S of the B s — > Di*" >± K T channels can be related to their counterparts x d and 5 d of 
the B d -> D ( * )=t 7r T decays as follows [2]: 



X d 



6, = 5a 



(27) 



These relations assume exact [/-spin symmetry; the impact of possible corrections will 
be addressed below. 

The BaBar |26| and Belle 27 collaborations have performed measurements which 
allow us to constrain the hadronic parameters \x d \ and 8 a . For the Bd —> D ± -n T system 
the following constraints have been extracted from studies of CP- violating effects [3] : 



,£>7I 



'lep 



-2\x d \ sin(0 d + 7) cos(5 d ) = -0.03 ± 0.017, 
-2\x d \ cos{(p d + 7) sin(<y = -0.022 ± 0.021. 



A corresponding analysis of the Bd — > D* ± ir T decays (for which L = 1) yields (3 



2|^|sin(0 d + 7) cos(<#) 



c^ = 2|^|cos(0 d + T ) sin(^) 



-0.039 ±0.010, 
-0.010 ±0.013, 



(28) 
(29) 

(30) 
(31) 



where we have used the label V to distinguish the vector D* system. In order to convert 
these experimental results into \x d \ and S d , we assume the value for 7 in (g with the 
B° d -B Q d mixing phase d = 2/3 = (42.8 ± 1.6)° g], which yields <\> d + 7 = (111 ± 7)°. 

Let us first extract \x d \ by determining the doubly Cabibbo-suppressed branching 
ratio BR(73° ->• D~n + ) from BR(73° 7J>77r+) with the help of the 517(3) flavour 
symmetry (28). Using the notation of Ref. [23], we write 



BR(B° 7J-7T+) = (±) BR(5° -> 7J77T- 



where 



N'pN'gN'w.. 



(32) 
(33) 



In analogy to (20), the $ are are phase-space factors, while 

r„ c-b> + / 2 n! 2 



[m 



D s 



f. 



D Fi 



and 



A/! 



Oi(D + 7T 

ai(7J+7r- 



(34) 



(35) 
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describe factorizable and non-factorizable S't/(3)-breaking effects, respectively. The M' E 
factor takes into account that B® — > D~tt + has a contribution from an exchange topology, 
which does not have a counterpart in the B® — > D~n + channel: 



ML 



T, 



D~ 



Td--k+ + E D - n + 
We then obtain the following additional constraint for Xd- 



(36) 




BR(5° -+ D-7T- 
BR(B° d D+7i- 



For the numerical analysis, we use the ratio of decay constants /d 8 //d 



(37) 



1.25 ± 



0.9771 ± 0.0009, where 



0.06 [24] and the form-factor ratio F 1 d (m 2 D )/F l 
we have applied the evolution equation for the B^ — > tt + form factor given in Ref. 129 



For the decays entering (32), factorization is not expected to work well. Indeed, follow- 
ing the approach discussed in Ref. 23 , we extract \ai(D+n~)\ = 0.68 ± 0.12 from the 
experimental data, while factorization would correspond to a value around one. Unfor- 



tunately, an analogous factorization test for B® — > D 7r + cannot be performe 



We 



allow for 20% 5'[/(3)-breaking effects for the non-factorizable contributions, i.e. for the 
deviation of |ai| from one, leading to M' a = 1.0 ± 0.2. 

In order to estimate the importance of the exchange contribution, we apply the SU(3) 
flavour symmetry and use experimental information on BR(5° — > D+K~) = (2.2±0.5) x 
10 -5 (24), which receives only contributions from exchange topologies. Comparing it to 
the contribution from tree topologies, which we fix again through BR(5^ — > D~ir + ) = 
(2.16 ± 0.26) x 10~ 5 (30|, we obtain: 



E 



T, 



D-tH 



fx 



ub 



V, 



cb 



/ BR( J B° D+K- 
BR(5S D 7 ^ 



0.1. 



(38) 



Consequently, we estimate M' E ~ 1.0 ±0.2. In comparison with the value of Me ~ 0.97 ± 



0.08 given after (23), this range is larger. Although the exchange topologies entering 



both quantities are estimated to have similar absolute size, the analysis performed in 



Ref. 23 indicates a large angle between the E and T amplitudes, which reduces the 

= (2.68±0.13)x 



impact of E on the amplitude ratio in Me- 

Using finally also the experimental branching ratio BR(Z?° — > D + n 
10 -3 24 , the relation in (37) gives 



\x d \ = 0.0163 ± 0.0011| BR ±0.0026| st/(3) 



0.0163 ±0.0028. 



(39) 



This value is consistent with the results for Xd given in Ref. 130] . Combining (39) with 
(28) and (29) allows, in principle, the determination of <frd ± 7 and S d up to discrete 



ambiguities. Unfortunately, a corresponding numerical fit leaves these parameters still 
largely unconstrained. 



4 The branching ratio quoted by the Particle Data Group 
would create a circular argument. 



24 



is constructed from (37), so using this 
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Figure 3: The confidence level contours for the \ 2 fit of the hadronic parameters \xd 
and Sd as discussed in the text, illustrating also the impact of the \xd\ constraint in (39) 



v\ 



5^ from the constraints in 
(|28h — (|3 1 h using a x 2 fit- For the former parameter set we also include the constraint in 



We proceed to extract the parameters and \x : 



(39). The fit gives the following results: 



\tJ - n oi66 +00025 

\^d\ — u - u±uu -0.0029' 



„Vl 



r+0.014 



\ x k I — 0.025_ OJHKS . 



("35^)° 
(146^)° 



(40) 
(41) 



where the errors give the 68% confidence level for each parameter. The xV n dof is 0.53 
and 0.00 for the non- vector and vector decays, respectively. In Fig. [3j we show the 
corresponding 39%, 68% and 95% confidence level regions in the f^H 3 ^! plane. Note 



that the constraint in (39) considerably reduces the uncertainty of the \xd\ parameter 
for the non-vector decay. 



Using (27), we hereby find 



X s - 0.311±o;o53|. nput ± °- 06 ls[/(3)' 



x^ - OA7toit\ inpnt ± 0.09| 5{/(3) , 



35 -4o|input ± 20 1 51/(3) 
146 -25| mput ±20| 5f/(3) 



(42) 
(43) 



where we allow for S'?7(3)-breaking effects of 20% for the x^ parameters and ±20° for 
the strong phases. In later applications of these results, the uncertainties associated with 
the x d V \ 5j parameters and the S'[/(3)-breaking effects will be combined in quadrature. 

Before using the hadronic parameters given above to predict the observables of the 
B s — > D^^K^ decays in Section [ij which serve as input for an experimental study, let 
us first discuss the extraction of <p s + 7 from these channels, with a special emphasis on 
multiple discrete ambiguities and their resolution. 
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4 Extraction of <f> s + 7 and Discrete Ambiguities 



For the extraction of S + 7 from the B s D s system, it is necessary to measure 

the following CP asymmetries from time-dependent, tagged analyses: 



a CP (B s (t) -> D^ + K' 



r{B° s {t) -> d { ; )+ k-) - r(B° a {t) -> 

C cos(AM s t) + 5 sin(AM s t) 
cosh(?/ s i/rsj + >W smh(y s tj r Bs ) ' 



(*)+&—> 



(44) 



an analogous expression holds for the CP-conjugate D^' K + final states, where C, S 
and A at are simply replaced by C, S and *4.Ar, respectively. The observables take the 
following form |2j: 



C 



1-xj 
1 + x\ 



C = + 



l + xl 



(45) 



S 



1) 



1 2X ~ sin(0 s + 7 + <U S=(-l) L -^sin(0 s + 7 -5 s ), (46) 



1 + x 



1 + x 



which complement the expressions for A&r and A&r i n (pf ■ 

For the following discussion, it is convenient to introduce the observable combinations 



as well as 



s+ = (-1) 



2 

2 Xe 



(5)+ = + cos 8 S sm((J) s + 7) 



1 + x 
2x s 



21 



sin8 s cos(0 s + 7), 



where 



X s 



I -(C)- 
1 + (C)- 



yielding 



1 + ^ 

2x x 



(47) 

(48) 
(49) 

(50) 



Finally, we obtain 



sin 



+ 7) = ±1 



(l + 4-S 2 )± ,/(l + 4- S 2 )2 _ 4s 2 



(51) 



which results in an eightfold solution for <p s + 7. This is the "conventional" extraction 
of this quantity [TJ|2] . 

As was pointed out in Ref. [2], the observable combinations 



(A A r)- 



Aaf + Aav 



2 x s 



cos 5,, cos(0 s + 7) 



(52) 
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(&+7)8 



(<t>B + 7)6 



Figure 4: Illustration of the complex numbers (Aat + iS) and (Aat + iS) with lengths 
vl — C 2 in the complex plane. Left panel: illustration of the conventional extraction of 
4> s + 1 and the associated eightfold discrete ambiguity (see (51)). Right panel: illustra- 
tion of the reduction of the discrete ambiguity to a twofold one through the untagged 
observables ^Ar and ^Ar (see (54)). 



(^Ar)- 



Aat — Aat 



2 x s 
l + xl 



sin 5 S sin(0 s + 7) 



(53) 



can be combined with the mixing- induced CP asymmetries (S)± to derive the relation 

(S)+ (A Ar )- 



tan( 



+ 7) 



(^Ar> + (S)- ' 

which allows the extraction of <p s + 7 up to a twofold ambiguity; moreover, we have 



(54) 



tan(0 s + 7)| 



(Aat) 2 - +(S)t 



'(A AT )l-(S): 



(55) 



(S) 2 - + (A Ar )l V (S)l -(A AT )Y 

The final ambiguity can be resolved from factorization arguments, where we expect 

cos5 s >0, cos5f<0, (56) 

a pattern that agrees well with the results of the £7-spin analysis presented in Section [SJ 
where the results for the strong phases in (42) and (43) give 



cos 5. 



0.821™ cos 5] 



-0.83lg.-g. 



(57) 



Combining this with (48), the sign of sin(0 s + 7) can then be determined. Thus, under 
reasonable assumptions, the extraction of S + 7 is unambiguous. 

It is instructive to illustrate these features, which can be hidden in a global experi- 
mental fit (see Section [5]). As the observables satisfy 



C 2 + S 2 + A 2 Ar 



1 = C 2 + S 2 + A 2 



(58) 
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and are hence not independent, only two of the three observables for each of the final 
states of the B s — > D i 's" >± K T system are needed for the determination of <p s + 7. We 
introduce the complex numbers^ 



A Ar + i S = -(-1) Vl - C 2 e -^+^) (59) 
A Ar + iS = -(-l) L \]l _^V^+7-<5 s ) ; (60) 



which, as C = —C (see (45)), have the same absolute value and thus span the same 



circle in the complex plane. The weak phase <p s + 7 corresponds to the polar angle of a 



complex number that lies exactly between (59) and (60), with an equal angular distance 



of S s to both. Let us first have a look at the conventional strategy, which does not use 
the information provided by the untagged A A r, A A r observables. The C and C then 
fix a circle in the complex plane, while the mixing-induced CP asymmetries S, S fix the 
component in the imaginary direction. As illustrated in the left panel of Fig. |4| this 
results in an eightfold discrete ambiguity for <j) s + 7. On the contrary, as shown in the 
right panel of Fig. [4], if the mixing-induced CP asymmetries are measured together with 
the untagged observables, the discrete ambiguity is reduced to a twofold one, which can 
be fully resolved as discussed above. Consequently, the optimal observable sets for the 
extraction of 4> s + 7 are S, S and A A r, 

Another important advantage of these observables is not only that they depend lin- 



early on x s - in contrast to C, C and the determination of this parameter through (25) 



- but that x s drops out in (54) and (55). Interestingly, as we will see in the next section, 
both observable sets can be accessed with similar precision at LHCb: the extraction of 
the untagged -4.Ar, A^r observables relies on the B s decay width parameter Q, while the 
measurement of the S, S observables requires the tagging of the flavour of the initially 
produced B® or B® mesons. 



5 Experimental Prospects 

The hadronic parameters determined in Section |3j with the phases in Q and ^ , allow 
us to make predictions of the observables of the B s — > DfK T decays: 

r eff = 0.971±85g t Bb , ^ = -0.491°;^ C = -0.8241°$?, 5 = 0.29^, 

7* = 1.025j£°» r Ba , A Ar = 0.1lt° o it C = 0.824t° °ll S=0.55j££. (61) 
In analogy, for the B s — > D* decays we obtain 

t£ = 0.954lffi r Ba , Al r = -0.66t il, C v = -0M%%, S V = 0.4S£™, 

fI S - 1-0271^ r Bs , A V Ar = 0.13±S^g, C V = 0M±°Z ^ = 0.76^. (62) 



5 Their relation to the complex observables £ and £ defined in Ref. 21 is given by 2£/(l + |^| 
Aaf + i S and 2£/(l + |^| 2 ) = ^ A r +iS, respectively. 



2\ 
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Table 1: Statistical uncertainties of B® — > DfK T CP observables for various data sam- 
ples as determined from our toy study. The difference in sensitivity of (*4.Ar) + , (^Ar)- 
is due to a correlation between A&r and ^.Ar of 0.5 observed in our toy simulations. 



Scenario 


a(C , (J) 


a(S , S) 


v{A A r , Aat) 




^((^Ar)+) 




LHCb 2012 


±0.176 


±0.252 


±0.210 


±0.173 


±0.194 


±0.113 


LHCb 2018 


±0.077 


±0.110 


±0.092 


±0.076 


±0.085 


±0.049 


LHCb Upgrade 


±0.032 


±0.046 


±0.038 


±0.032 


±0.035 


±0.020 



Table 2: Experimental uncertainties on the weak phase S + 7, strong phase 5 S and 
hadronic parameter x s for various data samples as determined from our toy simulations. 
Results for the "conventional method", which excludes the untagged observables, are 
also shown. The errors correspond to the central values (j) s + 7 = 65.5°, 5 S = —35° and 
x s = 0.31. 





With Air and ^.Ar 


"Conventional Method" 


Scenario 


4>s + 7 $s X s 


4> s + 7 S s x s 


LHCb 2012 
LHCb 2018 
LHCb Upgrade 


[±17]° [±17]° ±0.080 
[±7.3]° [±7.3]° ±0.035 
[±3.0]° [±3.0]° ±0.015 


±0.11 

[±%)° [±lt)° ±0.048 



Furthermore, our predictions for the branching ratio observables ([5]) and (13) are 

BR(B S -> DfKT) cxp 



BR(B S -> D+K~" 



/cxp 



BR(B S -> D;K+) exp 



BR(5 S -> DfK~ 



/cxp 



BR(5 S -> L>ji^ 



/cxp 



S£/(3) 



S£/(3) 



-0.027; 



-0.052 
-0.019) 



(63) 
(64) 



respectively. The prediction in (63) is compared to the current experimental results in 
Fig. [2j Similarly, we predict for the vector decays: 

BR(B S L> s *±i^) exp 



BR(B S -> 



/exp 



BR(5 S -> D s *±7rT) exp 
BR(5 S -> D*~K+) exp 



BR(5 S -> D s *+ii"-)exp + BR(B S -> Dr^ + ) 



cxp 



51/(3) 



51/(3) 



n n99 +a030 

u - uyy -0.036) 



-0 ns5 +0 - 056 



(65) 
(66) 



To estimate the experimental sensitivity for the observables, a simple Monte Carlo 
simulation has been performed, using as theoretical input the central values x s = 0.311, 

17.72 ps" 1 



31 , y s = 0.088 (see @) and 7 + 



65.5° 



S s = -35° (see (42)), Am, 
(see ([T]) and (J3|). A global fit to the decay distributions then simultaneously determines 
the observables given in (61). 



13 



l.OOr 
0.75' 
0.50 
0.25 - 
= 0.00- 
-0.25 
-0.50 
-0.75 

-l.oo L 



X-| 




[A 


M>+/ 




sol 
















¥#-8 + 7)«a 2 


— Central Value 
-- 68%C.L.(1D 
90%C.L.(1D 

— 39%C.L.(2D 
-- 68%C.L.(2D 





-0.75-0.5O-u. 25 u.l .25 0.50 0.1 

Tie 




-0.25 

— Central Value 

-0.50] 08%C.L.(1D) 

90%C.L.(1D) 

-0.75 — 39% C.L.(2D) 

-- 68%C.L.(2D) 

-1.00 



-0.75-0.50-0.25 0.01) 0.25 0.50 0.75 1.00 
Tie 



1.00 
0.75 
0.50 
0.25 
= 0.00 
-0.25 
-0.50 
-0.75 
-1.00 

1.00 
0.75 
0.50 
0.25 
= 0.00 
-0.25 
-0.50 
-0.75 
-1.00 




Central Value 

1% C.L.(ID) 
90%C.L.(1D) 
39%C.L.(2D) 

l%C.L.(2D) 

-0.50-0.25 0.00 0.25 0.50 0.1 
Tie 




Central Value 
i% C.L.(ID) 
90%C.L.(1D) 
39%C.L.(2D) 
68%C.L.(2D) 



-0.75-0.5O-D.25U.0ll 0.25 11.50 0.75 1.00 
Tie 



1.00 

0.75 

0.50 

0.25 

j 0.00 
1 

-0.25 
-0.50 
-0.75 
-1.00 

1.00 
0.75 
0.50 
0.25 
: 0.00 
-0.25 
-0.50 
-0.75 
-1.00 





[A 




((/>» + 7)«ol; 










V [4>, + 7)sol 2 


— Central Value 
-- 68%C.L.(1D 
90%C.L.(1D 

— 39%C.L.(2D 
-- 68%C.L.(2D 





-0.75-0.5im. 2", ii 25 0.51) 0.75 1.00 

Tie 



(A & r)- 



Central Value 
68% C.L.(ID) 
90% C.L.(ID) 
39% C.L.(2D) 
68% C.L.(2D) 




-0.75-0.50-0.25 0.00 0.25 0.50 0.75 1.00 
Tie 



Figure 5: Illustration of the determination of S + 7 from the separate observable com- 



binations (*4.Ar)+ and (S} + (top row) and (Aaf)- and (S)- (bottom row); see (67) and 



(68). The increasing experimental sensitivity of the panels from left to right corresponds 
to expectations of the LHCb experiment by the end of 2012, before the upgrade and 
after the upgrade, respectively 




<t>> + 7 [deg] </>., + 7 [deg] <t> a + 7 [deg] 



Figure 6: Illustration of the determination of S + 7 from (*4.Ar)+ and {S) + by means of 
(54). We also show the eightfold solution resulting from the "conventional" method (51) 
as discussed in the text. The increasing experimental sensitivity of the panels from left 
to right corresponds to expectations of the LHCb experiment by the end of 2012, before 
the upgrade and after the upgrade, respectively. 
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Figure 7: Illustration of the determination of S + 7 from a simultaneous fit to *4.Ar, S, 
C and their CP conjugates in the complex plane (see also Fig. [5]). The 68% confidence 
levels for (p s + 1 and the CP observables are indicated by the hatched and shaded regions, 
respectively. The two solutions shown for S + 7 and 8 S correspond to the remaining 
twofold discrete ambiguity discussed in the text. The increasing experimental sensitivity 
of the panels from left to right corresponds to expectations of the LHCb experiment by 
the end of 2012, before the upgrade and after the upgrade, respectively. 



Experimental data sets are simulated assuming approximate detector performance, 



as discussed in Ref. 31 by the LHCb collaboration, corresponding to a decay-time 
resolution of 50 fs, a flavour tagging efficiency of 38%, and a wrong-tag probability of 
34%. The sensitivity is estimated for data sets that would correspond to about 1100 
events per fb _1 of collected integrated luminosity [14] , selecting only B s candidates with 
a lifetime of t > 0.5 ps. Systematic effects, such as the presence of background events, 
are ignored in this study. 



In the toy simulation, the observables listed in (61) are determined from a fit to 
the decay distributions from 3500 simulated B s — > D^K T events corresponding to the 
approximate data sample that can be collected by the LHCb experiment by the end of 
2012. The fit is repeated for 2000 different data sets, resulting in an estimate for the 
sensitivity for the observables, which is comparable to the accuracy of the prediction 
itself. In Table [TJ the statistical uncertainties for the observables are listed for data 
samples corresponding to the expected integrated luminosity of the LHCb experiment 
at the end of 2012, before the upgrade, and after the upgrade. In our toy simulations an 
average correlation of 0.5 was observed between the *4.Ar and A&r observables, which is 
taken into account in the fits below; the correlations between the other CP observables 
is found to be negligible. 

As a final step, these estimated experimental uncertainties for the observables Aat, 
S, C and their CP conjugates can be translated into a determination for 7. Using only S, 
S and C, C following the "conventional" approach described by (51), the experimental 
sensitivity is not sufficient to determine 7 for a data sample of about 3500 events, which 
can be collected by the end of 2012. A factor five increase in data size, corresponding 
to the end of the current LHCb experiment, would result in a sensitivity of 7 + <p s = 
(65.6^26)° ^ the solution around the input value for 7 is selected. 
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If instead the observable pairs 



(A A ^) + + ^(S)- 



-(-V. 



1 - (C) 2 _ cos5 s e 



s+7) 



(67) 



and 



(A A r}-+i(S}- = -(-l) L y / l - (C) 2 _sm5 s e^ /2 -^ + ^ 



(68) 



are used separately, the 2012 data sample corresponds to experimental sensitivities for 
7 + 4> s of ±24° and ±29°, respectively; see the left panel of Fig. [5j In Fig. [6j we illustrate 
the extraction of 7 + (f) s from (S) + and (Aav)+ by means of the first relation in (54). 
Finally, combining all the observables, i.e. *4Ar, <S and C with their CP conjugates, 
improves the sensitivity to ±17°, as illustrated in the left panel of Fig. [7j 

With increasing data samples, shown in the middle and right panels of Figs. |5_ 
the precision on the measurement of 7 + S is expected to increase to about ±7° (3°) 
using 18k(130k) B s — > DfK^ events, which could finally be collected by the current 



(upgraded) LHCb experiment 32 , assuming unchanged trigger and tagging performance. 
The different experimental errors for the determination of <p s + 7, 5 S and x s from the 
B s — > DfK T decays are collected in Table [2j 

The magnitude of *4.Ar + iS can be further constrained through the SU(3) flavour 
symmetry, i.e. through (25) or by means of (27) with (37). However, we find that this 
input, which would introduce the SU(3) flavour symmetry into a theoretically clean 
strategy, does not significantly improve the precision for 7 + <ft s . 

On the other hand, if also decays of the type B s — > D*^K T can be reconstructed, 
the precision could be further enhanced in a theoretically clean way. These channels 
require the reconstruction of a radiative photon in the decay D* — > Df^ and as such 
are experimentally more challenging. In Ref. 33 , a gain in statistics of 28% is deemed 



possible, leading to an improvement of 13% on the determination of 4> s + 7. 



6 Conclusions 

The decays B s — >■ D s K T offer an interesting playground for the LHCb experiment in 
this decade. We have performed a detailed analysis of the observables of these channels, 
addressing in particular the impact of the sizable B s decay width difference Ar s , which 
has recently been established. This quantity leads to a subtle difference between the 
experimental and theoretical branching ratios of the B s — >■ D s K T decays, which can 
be resolved experimentally through time information on the corresponding untagged 
data samples, such as measurements of the effective decay lifetimes. We derived a lower 
bound for the ratio of the experimental B s — > DfK T and B s — > Dfir T branching ratios 
given in ([5]), and observe that the central value for the LHCb result is too small by about 
two standard deviations. 

The width difference AT S offers the untagged observables and ^4.Ar for the 

final states Dg + K~ and Dl*^~K + , respectively, which can nicely be combined with 
the corresponding mixing-induced CP asymmetries S and S to determine <fi s + 7 in 
an unambiguous way. We have illustrated this strategy and have obtained predictions 
for the B s — > Ds^K^ 1 observables from an SU(3) analysis of the I?-factory data for 
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Bd — > D^^it 11 , Bd — > Dfn T decays. Moreover, making experimental simulations, we 
have shown that the interplay between the untagged observables Aaf, ^Ar and the 
tagged CP asymmetries S, S is actually the key feature for being able to measure S + 7 
through the B s — > D^^K 11 decays at LHCb. In this sense, the favourably large value 
of Ar s is a present from Nature. 
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